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Abstract 

Background: Tliis paper describes tine "EMG Driven Force Estimator (EMGD-FE)", a 
Matlab® grapliical user interface (GUI) application that estinnates skeletal muscle 
forces from electromyography (EMG) signals. Muscle forces are obtained by 
numerically integrating a system of ordinary differential equations (ODEs) that 
simulates Hill-type muscle dynamics and that utilises EMG signals as input. In the 
current version, the GUI can estimate the forces of lower limb muscles executing 
isometric contractions. Muscles from other parts of the body can be tested as well, 
although no default values for model parameters are provided. To achieve accurate 
evaluations, EMG collection is performed simultaneously with torque measurement 
from a dynamometer. The computer application guides the user, step-by-step, to 
pre-process the raw EMG signals, create inputs for the muscle model, numerically 
integrate the ODEs and analyse the results. 

Results: An example of the application's functions is presented using the quadriceps 
femoris muscle. Individual muscle force estimations for the four components as well 
the knee isometric torque are shown. 

Conclusions: The proposed GUI can estimate individual muscle forces from EMG 
signals of skeletal muscles. The estimation accuracy depends on several factors, 
including signal collection and modelling hypothesis issues. 

Keywords: Muscle biomechanics, EMG-driven models. Muscle force estimation. 
Lower limb muscles 



Background 

Estimating skeletal muscle forces in vivo is a challenging problem in biomechanics. 
The possibilities for directly measuring such forces are limited due to the invasive na- 
ture of the procedure. Inverse multibody dynamics, associated with static or dynamic 
optimisation, can produce an estimate of a muscles dynamic state [1]; however, the re- 
sults are strongly dependent on the solution of the muscle force-sharing problem [2]. 
EMG-driven modelling is an attractive technique to estimate muscle forces in vivo for 
both healthy and pathological conditions [3-12]. The basic idea consists of using elec- 
tromyography (EMG) signals as the inputs of a Hill- type muscle model. A specific 
characteristic of the EMG-driven model used in EMGD-FE relies on that it is formu- 
lated entirely in terms of ordinary differential equations (ODEs) to represent muscle 
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dynamics [6]. Numerical integration of the ODEs yields the dynamic state of the muscle 
model In our case, it comprises the muscle active state (activation), the tendon force 
and the length of the contractile part of the musculotendon unit. 

Compared with force sharing solutions obtained through optimisation, EMG-driven 
models do not require the definition of an a priori cost function and they have a low 
numerical cost and converge easily, especially when compared with dynamic optimisa- 
tion [1]. In addition, EMG-driven models take into account muscle dynamics, which is 
usually neglected in static optimisation [13]. However, EMG techniques present some 
well-known limitations that apply to EMG -driven models as well. Deep muscles are not 
accessible by surface electrodes, and large-scale problems require numerous EMG 
channels. The relationship between the collected bipolar EMG and the real degree of 
muscle activity is not trivial, and it can be subjected to a number of misinterpretations 
[14]. Considering such limitations, bipolar EMG is still a very useful tool to assess the 
degree of muscle activity. In the best case, a rectified, MVC-normalised and filtered 
EMG gives a reasonable estimation of muscle excitation timing and amplitude, but it 
provides no force information. Considering muscle mechanics models further extend 
the utility of an EMG analysis, allowing the estimation of muscle force. 

Comprehensive predictive simulations — such as the solution of a dynamic optimisa- 
tion/optimal control problem— of muscle activity during a motor task cannot be per- 
formed by EMG -driven models; however, forces estimated using experimental EMG 
data can be used to constrain the solutions of the optimisation problem to a more 
physiologically feasible space. 

This paper presents a graphical user interface (GUI) that estimates muscle forces 
using a dynamic Hill-type EMG-driven model (EMG Driven Force Estimator (EMGD-FE) 
vl.O). It is implemented as an open source Matlab® graphical user interface (GUI) 
under general public license (GPL). Essentially, it automates the process that has been 
used by our group to process EMGs and to pass them into a dynamic model of skeletal 
muscles, which is then numerically integrated. The resulting muscle forces are multi- 
plied by their respective moment arms, and the partial joint torques are totalled and 
compared to a joint torque curve, obtained by inverse dynamic analysis or dynamome- 
try. Finally, the results are plotted and analysed^. All the processing steps are performed 
through sequential and self-explanatory GUI windows, which allow the next step to be 
performed only when the previous one is accomplished, thereby reducing mistakes. A 
set of default parameters for the lower limb is supplied and can be easily edited in a 
set of Extensible Markup Language (XML) files. Up to 6 Generic Muscles, without 
pre-defined parameters, can be used to analyse unlisted muscles. 

Several muscle-mechanics formulations can be used for EMG-driven studies. Here, 
we follow the analytical and the numerical procedures used in our own muscle bio- 
mechanics research, which include model formulation, parameter estimation, data col- 
lection and filtering characteristics. In the current GUI implementation, no parameter 
optimisation algorithm for the adjustment of torque curves is available. 

An example involving quadriceps femoris is included for illustrative purposes. Our 
intention is to provide periodic updates to the GUI that incorporate new features and im- 
provements into the muscle dynamics model, the EMG collection steps and the process- 
ing procedures. One major limitation of the current version is that only isometric analyses 
can be performed. For non-isometric muscle operation, the model is particularly sensitive 
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to certain parameters, such as the tendon slack length, the maximum force and the pen- 
nation angle. From preliminary simulations [15], we observed that in vivo calibration 
using both numerical optimisation as well as more sophisticated muscle and tendon ma- 
terial constitutive equations are necessary to obtain reliable force estimates. This feature 
will be included in future versions of the GUI. 

Implementation 

Muscle dynamics, data collection and EMG to force processing 

The Hill-type dynamic muscle model comprises contractile, parallel elastic and parallel 
damping elements, as well as an elastic tendon. It is a modified version of the Zajac 
[16] model, incorporating parallel elastic and damping elements with the main object- 
ive of improving numerical stability. It has three dynamic states, activation {a), tendon 
force (F^) and muscle length {L^), and two inputs excitation {u) and musculotendon 
velocity (V^^). The last input variable was set to zero because the contractions are iso- 
metric. The ODEs (Eq. 1) are integrated using a variable step-size Runge-Kutta algo- 
rithm. The first line of Eq. 1 represents the activation dynamics, and the second 
represents the contraction dynamics. The last line allows the integration of the muscle 
velocity to explicitly produce the muscle length. This state variable is used in a polyno- 
mial function that is part of the contraction dynamics, which relates the muscle max- 
imal force and length. The model has been used extensively in a number of published 
studies, and its analytical derivation is found in Menegaldo and Oliveira [7]. 



The EMG -driven force estimation procedure comprises the following: collecting 
EMG signals, followed by rectifying, filtering and numerically integrating the muscle 
model based on this input. The resulting tendon force is multiplied by the muscle mo- 
ment arm, and the sum of the individual muscles' contributions to the total torque can 
be compared to the simultaneously collected dynamometer measurements (Figure 1). 
Moment arm values depend on the joint angle that the isometric experiment was per- 
formed, and should be provided by the user. They can be easily estimated using Opensim 
[17] or polynomial regression equations [18]. 

In the current implementation, the activation dynamics comprises two steps: a differen- 
tial bilinear equation and an algebraic non-linearisation operation ("A-model", Figure 2), 
which is modulated by a single "A" parameter [19]. This model slightly augments low- 
level activations relative to high-level ones, representing relationships between experimen- 
tally observed EMG and isometric forces [20]. The effect on the prediction accuracy, by 
changing the "A" parameter in the EMG-driven model, for the quadriceps muscle has 
been addressed in Menegaldo and Oliveira [21,22]. 

The muscle contraction protocol, such as ramp, step, or sine, is freely chosen by the user 
provided that the contractions are isometric. Electrodes can be either bipolar or multi- 
channel, but a single raw EMG signal should be provided for each muscle. The maximum 
voluntary contraction (MVC) torque and EMG, as well as the data of the submaximal 
tests, should be collected without changing the electrodes or the subjects position. 



f(x,u) 




(1) 
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Figure 1 EMG-driven model for muscle force estimation. The EMG for muscles Mi 23 (or more) is 
rectified and low-pass filtered for envelope extraction. Neuromuscular excitation, u(t), is the input for the 
activation dynamics. Activation, a(t), drives the contraction dynamics, which is integrated to produce the 
tendon force. The joint torque, which is compared to the joint torque measured by the dynamometer, is 
estimated by the sum of the tendon forces multiplied by their respective moment arms (either positive or 
negative, depending on if the muscle actuates as an extensor or a flexor). 



Procedure of analysis using the EMGD-FE 

The procedure for using EMGD-FE is summarised in Figure 3 and is described below. The 
user is guided to provide the necessary information and to choose the appropriate analysis 
parameters for each case, such as the digital filter design, the selection of the signal time 
window that will be considered, and the muscle model parameters. Screenshots from se- 
lected windows are shown to illustrate the applications interface (Figures 4, 5, and 6). 
More details can be found in the EMGD-FE user s manual. 

1. Provide basic information about the subject and the collected data. A struct variable 
".emg" is created; 




Figure 2 Activation dynamics. Its implementation comprises two steps. The first is the second row of 
eq. (1), and the second is a non-linear algebraic equation, which augments low-level activations and is 
modulated by a single "A" parameter [19]. 

V ) 



Menegaldo et al. BioMedical Engineering Online 2014, 13:37 
http://www.bionnedical-engineering-online.conn/content/1 3/1/37 



Page 5 of 1 2 



1 .Create new 
subject /test 




r 

8. Non-measured 
muscles estimation 

L. 




9. Submaximal 
tests offset 
removal 


i 






2. Add MVC tests 




7. Select the 
muscles to be 
included in the 
analysis 




10. Set muscle 
model parameters 


_l 


_l 


_l 


3. MVC filtering 
and pre-processing 




6. Submaximal 
EMG testspre- 
processing 




11. Integrate 
muscle ODEs 


_i 


_l 


_l 


4. MVC offset 
removal 




5. MVC interval 
selection 




12. Analysis of the 
results 



Figure 3 Procedure implemented in tlie EMGD-FE GUI for estimating the forces based on 
EMG signals. 



2. A text file containing the MVC EMG and torque data must be provided. The data 
can be organised into a single or separate file(s). The MVC signal is expected to 
contain an initial timespan of a relaxed EMG (for a few seconds), an intermediate 
range with the maximum voluntary activation (approximately 5 seconds to avoid 
fatigue) and an additional relaxing timespan. Within this range, the torque signal 
has to be sustained for about 2 seconds. Curve shape with a sharp EMG peak 



I Set Torque Offset 



Drag the cursor over a flat region on the graph to calculate the offset of torque signal automaticaly 
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Figure 4 Selection of the maximum torque region to normalise EMG inputs by MVC EMG. 
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Figure 5 Estimation of thie unmeasured EMG signals of deep muscles through the application of a 
user-defined formula that may depend on the EMG signals from other muscles. 



should be avoided. These intervals are useful for setting the initial conditions of the 
ODEs' states and checking for the quality of the EMG baseline; 
3. Rectify and apply zero-lag (forward and backward) filters to the MVC EMG signal: 
band-pass to remove artefacts, mains hums filter^ (50 or 60 Hz) and its three first 
harmonics) and a low pass filter to extract the signal envelope. The joint torque sig- 
nal, measured through a dynamometer synchronised with the EMG, is also 
smoothed by the digital filtering. EMG and dynamometer sampling frequency 



I Model Parameters 

- Global Parameters 

Maximurr Tension [N/cm^] 
Activation Time [s] 
Deactivation Time (s] 
Norm. Tendon Stiffness x Lst 
PE Damping Constant [N.s/m] 
PE Stifness 

Contraction Velocity [m/s] 
Time Scale [s] 



I 30.3 
0.012 
0.024 
I 37.5 
I 225 
I 1500 



- Correction Factors — 

Max. Force 
Optimal Fiber Length 
Moment Arm [m] 



- Muscle Paramelers- 



Rectus Femoris 
Vastus Medians 
Vastus Intermedius 



[7] Set Maximum Force directly [N] 

Physiological Cross-Sectional Area [cm^] 

Moment Arm [m] 

Optimal Fiber Length [m] 

Optimal PennatiDn Angle [deg] 

Tendon Slack Length [m] 

Relative Initial Length of the Muscle 



I 1871 



I 0 048 



I Back I I Finish 



PI 
P2 
P3 
P4 
P5 
P6 
P7 
P8 
P9 



Figure 6 Muscle model parameter editing window. Default parameters defined in an XML file are 
uploaded automatically and can be changed by the user. In the right-side panel, generic parameters are 
defined and made available to the 'model.m' file. This function is integrated using the Runge-Kutta algorithm. 
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should be compatible with EMG hardware low-pass filter, i.e. greater than twice ana- 
log filter cut-off frequency. Several digital filter parameters can be setup by the user; 

4. The user visually selects a signal interval in which activation is approximately zero. 
Its average EMG level is subtracted from the entire signal, providing a DC-level 
compensation; 

5. The user selects another MVC reference signal region such that the average EMG 
level of this region is used to normalise the subsequent sub-maximal EMG inputs.. 
This epoch should last approximately 2 seconds, such that the torque value is stable 
in the vicinity of its maximum plateau [23], as shown in Figure 4; 

6. The MVC procedure is repeated for the sub-maximal tests. The rectified and 
filtered sub-maximal EMG signals are normalised by the processed MVC EMG, as 
described above; 

7. Select the muscles, which cross the same joint, to be included in the analysis; 

8. Estimate deep muscles with unmeasured EMGs by providing a user-defined formula 
that may depend on the EMG signals from other muscles (Figure 5); 

9. Remove the offset. The user selects a signal range with zero activation, which will 
be subtracted from the whole signal. The user must be sure to record a portion of 
the EMG (approximately 5 seconds) with the muscles relaxed; 

10. Set muscle model parameters. Default values, using data from the Opensim lower 
limb model [17], are loaded from a XML file. A window shows the nominal 
parameters for each muscle that the user is allowed to change (Figure 6). Scaling 
factors can be set for certain parameters (left lower part of the window). In the 
right-side panel, a list of arbitrary parameters (P1-P9) can be defined, which will be 
passed to the Matlab m-function where the muscle model is defined. It allows 
some flexibility for the user changing the model equations without altering the 
GUI structure. In the current distribution, the parameters PI, P2 and P3 are used 
to set fixed or variable pennation angle relationships with joint angle, only for 
triceps surae muscle (see EMG-FE User Guide, Sec, 2.6.1. Setting Muscle 
Parameters, for details). The other parameters (P4-P9) have no use in the current 
version. Additionally, the activation dynamics "A" parameter is defined; 

11. Integration of the muscle dynamics ODEs. The model is written in a Matlab 
m-function called by the "ODE-45" Runge-Kutta integrator. The function where the 
model is defined has simple syntax, and the user can easily define or modify his/her 
model formulation. The memory requirements and the simulation time are small 
based on the fact that the simulations run in only a few seconds on an ordinary lap- 
top. Alternative muscle dynamics formulations that lead to stiff equations can take 
advantage of an appropriate ODE integration algorithm. Matlab provides a number of 
alternative integrators, which can easily replace the ODE45 using the same syntax. 

12. Several output options are available: the raw and processed EMG signals, the 
muscle forces, the muscle contributions to total torque, and the torque 
reproduction error. 

Results and discussion 

Sample results from an analysis using the quadriceps femoris muscle have been in- 
cluded for illustrative purposes. Data from vastus medialis (VM), vastus lateralis (VL) 
and rectus femoris (RF) were acquired at a 2 kHz sampling frequency in a young adult 
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male. Figure 7 shows the experimental setup used to obtain the bipolar EMG signals 
and the isometric knee torque. 

The isometric contractions were performed with the knee flexed at 90° (see details in 
[21,22]). Vastus intermedius (VI) EMG activity was estimated as the average between 
vastus medialis (VM) and lateralis (VL), as in Figure 5. Figure 8 shows the individual 
contribution of each quadriceps component to the total knee extension torque in a 
step-following manoeuver: the subject tries to track a protocol mask that corresponds 
to a 30-second plateau of 20% of the MVC torque, preceded and followed by relaxing 
intervals. In Figure 9, the total knee torque estimated from the EMG signals is com- 
pared with the torque measured from the dynamometer. 

The accuracy of individual muscle force estimation from EMG-driven models cannot 
be directly assessed. The joint torque error between the dynamometer measurements 
and the sum of the estimated partial muscle-force contributions is the best available 
index. For this particular experiment, the average RMS error between the estimated 
and the measured torque was 4 Nm, or approximately 10% of the plateau value. 

The user should be aware that model estimation accuracy strongly depends on the 
choice of muscle parameters. The most influential is maximum muscle force, which 
varies widely among individuals. Some techniques can be applied to estimate this par- 
ameter individually (see details in [7,21,22,24]). The influence of 'A' parameter on the 
knee torque estimation accuracy was addressed in [21,22]. Optimisation of certain par- 
ticular muscle model parameters, such as the maximum force, fibre length and tendon 
length, is a useful tool and easily decreases the torque error; however, it does not en- 
sure that the correct' parameters values are always identified [15]. 

An important question regarding EMG-driven models is the estimation of a deep 
muscle from neighbouring superficial muscles, such as the vastus intermedius case 
shown here. Watanabe and Akima [25] studied the relationship among EMG activities 
from knee torque, vastus intermedius and other components of quadriceps femoris. 
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Figure 8 Sample test data from quadriceps. The torque contribution of eacli component is summed 
and compared to tlie total estimated l<nee torque. 



Based on their data, we adjusted regression equations to estimate VI activity from the 
normaUsed torque; from the VM-, VL- and RF-measured EMGs; and as the average be- 
tween VL and VM [22]. These regression equations were used to estimate the quadri- 
ceps components' forces using the EMG-driven model. All estimation methods 
provided the same level of torque reconstruction accuracy, suggesting that the average 




-10' ' ' ' ' ' ' ' ' ' ' 

0 5 10 15 20 25 30 35 40 45 50 

time(s) 

Figure 9 Sample test data from quadriceps - torque. Data are presented as a comparison between tlie 
EMG-driven-estimated and tlie dynamometer-measured torque. 
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between VL and VM is a reasonable estimation of VI EMG, at least for this particular 
application. Because the GUI can be applied to virtually any joint, the user must be 
aware of the hypotheses that are being taken into account to estimate the unmeasured 
EMG from synergistic muscles. 

Conclusions 

This paper describes a Matlab application with a GUI that guides the user, step-by- 
step, to estimate the muscle forces from EMG signals. The forces are found by integrat- 
ing a set of ODEs describing a Hill-type muscle model. The EMG envelope is consid- 
ered the muscle excitation input. The models accuracy is assessed by comparing the 
estimated joint torque, which is calculated from the sum of the partial torque contribu- 
tions provided by the muscles, with that measured using a dynamometer. The GUI also 
provides several types of outputs and analysis tools. Although lower limb muscle refer- 
ence values are provided, the choice of muscle dynamics parameters depends on user 
discretionary. Force estimation accuracy is greatly sensitive on muscle parameter values 
and EMGD-FE can be useful for tuning these parameters. 

A new version of the GUI is being developed to broaden the possible types of EMG- 
driven analysis, which are currently restricted to isometric tests. Non-isometric analysis 
requires describing the musculotendon length and the moment arms as a function of 
the joint angles. In addition, the Hill-type muscle model becomes especially sensitive to 
certain parameters, such as the tendon stiffness and the slack length, requiring previous 
calibrations based on optimisation methods. 

Availability and requirements 

Project name: EMG Driven Force Estimator (EMGD-FE) 
Project home page: to be defined 
Operating system(s): Windows 
Programming language: Matlab 

Other requirements: Matlab Signal Processing Toolbox 
Any restrictions to use by non-academics: no restrictions 

Endnotes 

^The EMGD-FE v. 1.0 package, including synchronized EMG and torque samples for 
quadriceps femoris and triceps surae and the User Manual can be freely downloaded 
from http://www.peb.ufrj.br/docentes/Luciano/EMG-FE.htm. Sample data for other 
muscles can be made available on a collaborative basis. 

'The Matlab routine for filtering main hums noise has been developed by the Laboratorio 
di Ingegneria del Sistema Neuromuscolare e della Riabilitazione Motoria, Politecnico di 
Torino, Italy and is used with permission. 

Abbreviations 

EMG: Electromyography; RF: rectus femoris; VM: vastus medialis; VL: vastus lateralis; VI: vastus intermedius; GUI: Graphical 
user interface; ODE: Ordinary differential equation; XML: Extensible Markup Language; MVC: Maximum voluntary 
contraction. 
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